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Introduction 


This report is ihe final report for Cooperative Agreement NCC 2- 
1105: “In Situ Measurements of N20 and CH4 during SOLVE on the ER-2 
using a new tunable diode laser instrument.” The tasks outlined in the 
proposal are listed below with a brief comment. The publications and the 
conference presentations are listed. Finally the important publications are 
attached. 

The Cooperative Agreement made possible a research effort to 
produce high precision and high accuracy in situ measurements of methane 
and nitrous oxide on the ER-2 during the SOLVE field campaign and to 
analyze these measurements. These measurements of CH4 and N20 were of 
utmost importance to studies of the ozone losses in the Arctic winter and 
spring. The concentrations measured over a large spatial and temporal range 
allowed the separation of the dynamical and chemical ozone loss. 

The most important results of the SOLVE program were contained in 
two scientific papers (attached). This Cooperative Agreement allowed the 
participation of the Argus instrument in the program and the analysis of the 
data. 

II. NCC 2 1105 Tasks and Accomplishments 

The main tasks relating to the participation of the Argus instrument in 
SOLVE were as follows: 

Task 1: Perform intercomparisons of the different instruments on 

balloon, DC-8, and ER-2 platforms measuring the same species 
as our instrument. 

Task 2: Assess the subsidence in the Arctic vortex, specifically in 

conjunc tion with measurements aboard the ER-2 payload. 

Task 3: Study the dynamical history of the air masses sampled during 

the mission, especially with theory groups who calculate 
trajectories. 


1 



Task 4: Investigate small scale structures observed in profiles to 

improve our understanding of exchange of intra- and extra 
polar vortex air. Previous studies focused mainly on the 
analysis of ozone depleted air observed outside the vortex. 

Task 5: Investigate seasonal, interannual and interhemispheric change 

and differences in N20 and CH4, using the long time series of 
ATLAS N20 measurements. 

Task 6: Study the N20:N0y correlations to investigate denitrification. 

Task 7: Quantify the 03 loss due to chemistry by helping to estimate 

the initial ozone using correlations with long live tracers (e.g., 
N20). 

Task 8: Use the CH4 measurements and H20 to investigate the 

occurence of dehydration in the Arctic. 


In addition the following necessary tasks relating to the maintenance of the 
instrument and the laboratory were performed: 

Task 9. Maintenance and upgrading of the Argus airborne instrument 

9.1 Maintain Argus calibrations (ISO 9000 protocols) 

9.2 Maintain Argus spares 

9.3 Maintain Argus field support infrastructure 

9.4 Define Argus upgrades 

9.5 Argus documentation (ISO 9000) 

9.6 GSE and software operational maintenance and upgrades 

Task 10. Specific tasks to support lab infrastructure 

10.1 Maintain ESD capability /training (ISO 9000) 

10.2 Maintain electronic schematics (OrCAD software) 

10.3 Integrate LabView data acquisition into lab laser operations 
optical bench 

10.4 Maintain lab safety plan/haz-mat handling (ISO 9000 and 
HazMat regulations) 
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III. Publications: 


1. D.W. Fahey, R.S. Gao, K.S. Carslaw, J. Kettleborough, P.J. Popp, 

M.J. Northway, J.C. Holecek, S.C. Ciciora, RJ. McLaughlin, T.L. 
Thompson, R.H. Winkler, D.G. Baumgardner, B. Gandrud, P.O. 
Wennberg, S. Dhaniyala, K. McKinney, Th. Peter, R.J. Salawitch, 

T.P. Bui, J.W. Elkins, C.R. Webster, E.L. Atlas, H. Jost, J.C. Wilson, 

R.L. Herman, A. Kleinbohl, and M. vonKonig, The Detection of 
Large HN03-Containing Particles in the Winter Arctic Stratosphere, 
Science, Feb 9 2001: 1026-1031. 

2. Popp, P.J.; Northway, M.J.; Holecek, J.C.; Gao, R.S.; Fahey, D.W.; 
Elkins, J.W.; Hurst, D.F.; Romashkin, P.A.; Toon, G.C.; Sen, B.; 
Schauffler, S.M.; Salawitch, R.J.; Webster, C.R.; Herman, R.L.; Jost, 

H.; Bui, T.P.; Newman, P.A.; Lait, L.R. ’’Severe and Extensive 
Denitrification in the 1999-2000 Arctic Winter Stratosphere”, 

Geophys. Res. Lett Vol. 28, No. 15, p. 2875,2001 

3. Jost, H., M. Loewenstein, Argus: A Lightweight TDL Instrument to 
Measure Stratospheric Tracers, Application of Tunable Diode and 
Other Infrared Sources for Atmospheric Studies and Industrial Processing 
Monitoring II, SP1E annual meeting, Denver 1999. 


IV Presentations 

1. H. Jost, et al, “Mixing events during SOLVE/THESEO revealed by 
Tracer correlations” Annual Meeting of the Amer. Geophys. Union, 
December 2000, San Francisco, CA 

2. Jost, H., et al. , Evidence for transport events into the tropical middle 
Stratosphere, AGU spring meeting, Washington DC, 2000. 

3. Greenblatt, B.J., II. Jost et al,. Vortex Conditions in the Arctic Winter 
Stratosphere in 1999-2000: ER-2 Argus and Balloon LACE Tracer”, 
AGU spring meeting, Washington DC, 2000. 
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ABSTRACT 

Argus is a two channel, tunable diode laser instrument which measures atmospheric methane and nitrous oxide in the upper 
troposphere and stratosphere up to an altitude of 32 km using second harmonic detection. Investigations of stratospheric 
transport from mid-latitudes into the tropics is the current focus of our work which requires high precision and high accuracy 
data. Argus was designed for use on remotely piloted aircraft or light-payload balloons and weighs less than 20 kg. The two 
channels each have their own laser, optics, detector and signal processing chains. We sample methane at 3.3 urn and nitrous 
oxide at 4.5 at a rate of 0.1 Hz. The gas is sampled in a Herriott cell which has a total path length of 1 8.8 m. Each laser 
is current- and temperature-controlled by a dedicated microprocessor. We sweep the laser at 10 Hz and record direct 
absorption spectra during the ramping of the laser current. The lasers arc modulated with a 40 kHz sine wave, a phase 
sensitive amplifier and integrator detects the second harmonic data. The analysis is performed offline by applying direct fits to 
the measured spectra using the non-linear Marquardt-Levenberg algorithm. We' have recently flown Argus and ATLAS 
together on the ER-2 platform for comparison. Precision of Argus is 0.6% and accuracy is estimated to 4- 1 7% increasing with 
altitude. 

Keywords: atmospheric tracers, second harmonic detection, stratosphere, methane, nitrous oxide 


1. INTRODUCTION 

High accuracy and precision measurements of trace gases in the atmosphere allow us to investigate transport of air masses. In 
the stratosphere (10-45km above ground), relatively long lived tracers, like nitrous oxide (NjO), methane (CH 4 ), Carbon 
Dioxide (C0 3 ) or certain CFCs, contain information about the transport processes, because their lifetime is long compared to 
typical transport time scales. Both methane and nitrous oxide have sources in the troposphere and sinks in the stratosphere. 
At these altitudes the primary source of NjO is from air entering the stratosphere with a mixing ratio of about 3 15 ppbv with 
no annual cycle and a long term increase of 0.3% per year. In the stratosphere N 2 0 on average decreases monotonically with 
increasing altitude due to UV photolysis (wavelengths below 230 run) and because of its low chemical reactivity NjO has no 
other significant sink in the lower stratosphere. The lower stratospheric lifetime of N 2 0 is long compared to vertical and 
horizontal transport time scales 1 . Methane also has its primary source in the troposphere and enters the stratosphere with a 
mixing ratio of about 1.7 ppmv. Its decrease with increasing altitude is dominated by oxidation with OH and 0{ l D) with 
minor contribution from reaction with Cl. Its photochemical lifetime ranges from over a hundred years at 20 km to a few 
months at 40 km, which is long compared to the transport time constants (= weeks) 1 . Instead of investigating profiles of 
tracers, their correlation contains very useful information. The correlation of N 3 0 and CH 4 with other species is often found to 
be very compact or even linear. As has been explained by Plumb and Ko \ compact and linear relationships between tracers 
is found as long as the tracers are long lived relative to vertical transport time. The compact relationships last, although not 
linear, as long as the chemical lifetimes are still long compared to horizontal transport times. Inversely, if the correlation in 
the tropics differs from that measured at mid-latitude, the assumption of rapid mixing of the air masses seem to be violated 
and we can infer the effectiveness of the dynamical barrier, the so called tropical pipe 2,1 . Observation of laminae with 
exceptional correlations in the tropics can give us insight into the short term mixing of air masses from mid- latitudes into the 
tropics'. AH these investigations require high accuracy and precision measurements of trace gases. 

In search of a new high altitude research platform, the National Aeronautics and Space Administration (NASA) initiated work 
in the early 1 990 ’ s on the Perseus-A remotely piloted aircraft (RPA). As this platform had quite severe restrictions on 
payload weight and size, the need for small, lightweight in-situ instruments for atmospheric research came up. Argus, named 
in reference to the Greek multi-eyed creature because it measures several species (methane and nitrous oxide), was designed 
and realized for the Perseus-A vehicle Due to delays in the RPA project, Argus was deployed on a balloon payload. 
Together with five other tracer instrume nts it forms the Observation of Middle Stratosphere (OMS) payload. Within about 
Correspondence: Email: hjost@mail.arc nasa.gov; phone 650 604 0697; fax 650 604 3625 
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1.5 hours the 1 13,000 m 5 helium filled balloon rises to 32 km and is then slowly valved down before the payload is 
separated from the balloon and lands on a parachute. The total payload weighs about 600 kg and has been retrieved with 
minor or no damage after each of the seven flights. OMS has been deployed successfully from Northeastern Brazil, New 
Mexico and Alaska. 

Argus is an evolution from the very successful Airborne Tunable Laser Absorption Spectrometer (ATLAS) 5 . ATLAS is a 
single channel, diode laser instrument for nitrous oxide operating on the NASA ER-2 aircraft during many airborne field 
campaigns like the Airborne Antarctic Ozone Expedition (AAOE, 1987), or more recently the Stratospheric Tracers of 
Atmospheric Transport (STRAT, 1995*1997) and the Photochemistry of Ozone Loss in the Arctic Region in Summer 
(POLARIS, 1997). ATLAS measures the N 2 0 amount by forming the second harmonic peak ratio to a NIST traceable 
standard gas, whereas Argus records full spectra and uses spectroscopic line parameters to determine the gas amount. 

2. ARGUS INSTRUMENT DESCRIPTION 

2.1 Mechanical Configuration 

Argus in its current configuration, weighs about 20 kg and the main dimensions are 40cm x 30cm x 30 cm. The shroud, 
made out of fiber glass honeycomb and fiber glass facings, encloses the optics, electronics and sample cell. The pump, pump 
controller, calibration system and the platform interface box are mounted separately. A schematic is shown in Figure 1. The 
turbo pump is controlled to draw a constant flow of ambient air through the sample cell. The platform interface collects the 
Global Positioning System (GPS) signal provided by the balloon platform and transfers instrument data to the telemetry 
transmitter. 


turbo 

air inlet pump 



laser beam Argus box 

Figure 1: Schematic of Argus 

Inside the shroud, there are two main platforms. The lower level accommodates the electronics, electrical distribution, inlet 
gas heater and the MKS pressure transducer. Mounted on shock absorbers on top is the optical table, a honey comb 
composite structure. The aluminum dewar extends into the lower level and holds 1.8 liters of liquid nitrogen lasting for 10 
hours of operation and 30 hours idle. To reduce waiping of the table a custom shaped heater is mounted to the underside and 
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a PID control loop keeps the temperature at 30 °C. The electronics are kept at ambient pressure, as well as all the other 
components contained in the instrument box. Copper heat sinks have been glued onto the processors and other hot spots to 
prevent overheating due to the lack of convective air cooling at altitude where the air density is reduced to about 1% of an 
atmosphere. These heat sinks transfer the heat to the electronic chassis bottom and from there a 2.5cm aluminum bar connects 
to a radiation panel that also serves as a physical protection at landing. 

2.2 Optical Design 

The two channels each have their own laser, optics and detectors. Both channels have the same layout and are separated 
vertically. Only one is shown in Figure 1; they lie on top of each other in this perspective. We currently use lead salt diode 
lasers to measure methane in the v r band at 3.3pm and nitrous oxide in the v,-band at 4.5pm. The diodes, along with 4 
indium antimonide (InSb) detectors, are mounted in a single custom liquid nitrogen dewar. 

The f71 diode laser beam is shaped to a quasi-collimated £740 beam using two 1/2 inch diameter zinc sulfide lenses. At the 
intermediate focus is a pin hole to stop additional spatial modes a laser might have. All lenses and windows are anti- 
reflection coated with reflection less than 0.3%. 

A 10:1 beam splitter couples part of ihe laser beam through a reference cell to the dedicated reference detector. We use the 
well known spacing of absorption lines for the frequency calibration of the atmospheric spectra. The nitrous oxide channel 
reference cell contains CO and N 2 0 and we employ the CO R17 and N 2 0 P19 spacing which is 0.305 cm' 1 , the methane 
reference cell contains ,2 CH< and isotopic lJ CH 4 and we use the R0 spacing which is 0.100 cm respectively. Additionally, 
the reference spectral line shapes are used for instrument performance monitoring in flight. 

The beam passes 72 times through a 26 cm base path length Herriott cell for a total path of 18.8m. The ambient inlet air is 
preheated to 30*C. 

Currently the signal to noise of the instrument is limited by fringe patterns. The fringe spacing indicates reflecting surfaces 
separated by twice the Herriott cell length and we successfully reduced the fringe size by thermally modulating the walls of 
the Herriott cell, even though we were not yet able to identify the scattering sources that produce the dominant fringe. 

23 Electrical Design 

Figure 2 gives a simplified overview of the electronic design. We use the rapid scan method and second harmonic detection*. 
A first digital to analog converter (DAC) sets the bias current of the laser diode. We sweep the laser at 10 Hz and record 
direct absorption spectra during the ramping of the laser current. A ramp consists of 1024 discrete points controlled by the 
ramp DAC. To determine the detector dark current, the laser is chopped by turning the bias current off during the first 24 
steps of the ramp. A synthesized 40 kHz sine wave modulation is superimposed on the ramp by a third DAC. The sum of 
the 3 DAC signals drives a transistor that provides the laser current. The detector signal is amplified with a preamplifier and 
a programmable gain amplifier adjusts the signal within the dynamic range of the analog to digital converter (ADC). The 
phase sensitive amplifier demodulates at 80 kHz using a digital signal provided by the field programmable gate array 
(FPGA). An integrator is reset at the beginning of each ramp step and integrates during the four full cycles of the fundamental 
modulation frequency. The integrator signal is then digitized by an ADC. Due to the weight and size constraints we use 
stacked custom PCI 04 cards. Each laser is independently controlled by the FPGA, which initiates bias, ramp and 
modulation. Temperature is controlled by a channel-dedicated 386 microprocessor, which also coadds the spectra and 
communicates with the host processor. We attain an RMS temperature stability of about 5 mK using diodes as sensing 
elements. To maximize the 2f peak amplitude the sample modulation amplitude is continuously adjusted to 2.2 times the 
Voigt width of the absorption feature We coadd 100 signal ramps before storing them, corresponding to a 0.1 Hz sample 
rate. A full cycle of spectra consists of 2s coadded direct sample absorption, 10 s second harmonic sample absorption and 2s 
each of reference direct and second harmonic spectra. 

A third microprocessor, the host, running DOS operating system, controls the overall system and optics temperatures, serves 
an ethemet port which connects to the ground support computer, and saves the recorded spectra and housekeeping data to a 
solid state disk. Pressure in the sample cell is measured with two transducers: One with low accuracy for the range of 0-1000 
mbar, the other for 0-130 mbar with an accuracy of 1% of full scale. During the flight, the data is transmitted to the ground 
using a 19200 baud radio link. 

The electronics are exposed to ambient pressure and the heat is dissipated from a radiation panel connected to the electronics 
through a heat-conducting bar, as mentioned earlier. 
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Figure 2: simplif cd schematic of laser drive electronic and signal detection for one channel. 

2.4 Inflight Calibration system 

We have recently added an inflight calibration system to Argus that is based on a 1 liter gas tank Filled with a NOAA 
Climate Monitoring and Dynamics Laboratory (CMDL) certified whole air containing 313.6 ppb N 2 0 and 1766.3 ppb CH«. 
A regulator and a needle valve serving as a restrictor supply a flow that is twice as large as the flow provided by the pump. 
When the calibration system is activated by opening a solenoid, the cal gas flow back flushes the sample cell within a few 
seconds.This design makes sure that we don't get any contamination of the sample flow by leaks out of the calibration 
system, except if the solenoid fails. In this case the limited calibration gas supply will stop after a few minutes and 
undisturbed ambient measurements resume. The calibration system is activated 4 to 6 times for about two minutes during a 
flight by sending a start and a stop command over the telemetry. In the data analysis, the calibrations are not specially 
marked because the activation of the calibration system is independent of the instrument operation. We generally do not use 
the in-flight calibration data to correct the ambient concentrations, but to validate the full set of laboratory calibrations with a 
partial set during the flight and to estimate our accuracy. 

3. Instrument Calibration and Preparation 

In preparation for and after a flight we perform a calibration of the instrument at several pressures and mixing ratios. Whole air 
certified to 2% by NOAA/CMDL is flow mixed with blow off from a liquid nitrogen tank which proved to be a good zero gas 
for NjO and CH 4 . The flow mixing is controlled by measuring the flow of the calibration tank and using a metering valve on 
the zero gas to attain a certain mixirg ratio. A full set of calibration consists of 1 minute each of 100 %, 75%, 50%, 25% of 
whole air at a pressure of 1 30 mbar, 1 00 mbar, 70 mbar, 40 mbar, 25 mbar and 1 2 mbar. At each pressure we perform a linear 
least square fit to the retrieved volume mixing ratios to determine the calibration factor. Typical 1 a standard deviation of the 
slope is 1%. From the reference cel! data we determine the instrument contribution to the Gaussian broadening of the lines 
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f:l Ca rr Pare . U t0 5“ Iek ' CnCe ceil D °PP>er width retrieved during a flight. A saturation cell is inserted in the laser beam 

° f ? e 'T r and ' h l CUrrent and tempCra,UrC SC ‘ P 0 '" 15 ^ ad ^ USted f0r minima, 
h,n*Vk . ■ detune the temperature by as much as 0.3 K to check for laser mode breaks close to our desired lines and 

aser an7LT nl°heV entered K betWe ? With 3 3 Cm Iong 8 enTlanium cta,on we monitor the frequency tuning of the 

lsin P thi Zr SI r ear ;° better than ' percent > which allows us 10 derive the frequency axes of our spectra by simply 
using the spacing of two line absorptions in the reference cell. F y P y 

4. Data Analysis 

com z::7: r :t^ disk from r the instnj ™ nt r anaiyzc ,he da,a off - iinc - ^ ^ ^ i 00 ks an d 

the co iTes pond m eb free ^ s nertr T l catastroph ' c land,n 8' ln a ‘> rst ste P 'he 2f spectra are gain normalized and divided by 
me corresponding direct spectra. For frequency calibration we use spectra from the reference cells 

Jr' rC n C ‘ I s oflh e second harmonic component of a sine modulated Voigt function are applied to the preprocessed spectra usine 
the non-l inear Marquardt-Levenberg algorithm’. This method is smoothly vaiying between an inverse-Hessian (close to thf 
mmmium) and steepest descent (far from the minimum). We fit a total of 7 parameters: the .me porifio^ sTope and off*, of 

L ° rer l!f 3nd Dop P Ier width and ,he modulation amplitude. A typical specmirn Tnd fit is shown 
L oZler width 7 C °“ ld reducc thc sensitivit y residua. fringe patterns in Jspectra by fixing both Loren^ 

the HITRAN 96 snertr ca CU al * Lorentz w,dth fronl the measured pressure using the pressure broadening parameter from 
4^™, SpeC ^ OSCOp,C databa5e - For 'he Doppler width we use the mean value found by fitting the reference second 

fringe ratio ^rsMct^ with rhT!M CSSm r 8 ^ W “ h Z ^ Ch -‘ SqUare ’ 100 h '& number iterations, too small signal to 
The^algorithm is^mplemented in \ t ®. mp " ature change relative to the preceding spectrum are omitted from the final data, 
algorithm ,s implemented in Interactive Data Language (IDL™) and run on a Silicon Graphics workstation. 
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Figure 3: Top: measured (during the 970214 flight in Juaze.ro do Norte, Brazil. 7*8 and result of non-linear 

fit. Bottom. Difference between calculated and measured spectrum shown in top portion. 
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Figure 4: Profiles measured with Argus in Juazeiro do norte, Brazil 971111 (7*S), Fort Sumner, New Mexico, 980518 (31°N) and 

Fairbanks, Alaska, 970630 (64*N) 

5. Results 

Argus has been deployed on the OMS payload out of Fort Sumner, NM, on 960610, 960921 and 980518 (format 
YYMMDD), Juazeiro do norte, Brazil, on 970214, 971 111 and 971120, and Fairbanks, AK, on 970630. Data representing 
tropical data from Brazil, mid-latitude from Fort Sumner and high latitudes from Alaska are shown in Figure 4. The profiles 
are not very smooth, but have thin laminae (thickness less than 1 km) with depleted N 2 0 and CH 4 . In the Arctic, these 
laminae have been attributed to remnants of the polar vortex which broke up more than 2 month earlier 1 . In the tropical data, 
the laminae are recent intrusions of mid-latitude air and their age can be estimated from the ozone recovery to typical tropical 
values at the observed altitudes 4 . In Table 1 an overview of the error budget is shown for 9805 18. The largest contribution is 
from the pressure transducer (accuracy: 1 % of full scale ( 1 30mbar)) and can be up to 16% at 8 mbar, but of course decreases 



channel 1 

channel 2 


relative error 

relative error in 

Source 

in VMR (%) 

VMR (%) 

Gamma d 

1-3 

n/a 

Calibration Standard 

2 

2 

Calibration Slope 

i 

1,2 

4.2 

Pressure 

0.4-16 

0.4-16 

Total (RMS) 

2.6-16 

4.7-17 


Table I: accuracy error budget of Argus for 980518. 
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rapidly with increasing pressure. The calibration slope is determined from the pre- and post-flight laboratory calibration runs. 
For channel 2 we had a discrepancy in the two calibrations slopes and this reflected in the higher percentage error. Gamma_d 
reflects the uncertainty in the Gaussian line width contribution from the instrument. In November 98 we flew ATLAS and 
Argus together on the ER-2. Argus was run in a slightly different mode than on the balloon: For the signal second harmonic 
we co-added only 20 ramps before storing them, hence increasing the data rate to 0.5 Hz. The in-flight calibration system was 
not easily adaptable for the ER-2 and hence not used. 

From the preliminary data analysis, we found that Argus currently has a 1% precision in flight at 50 mbar for 2s data and 
0.6% for 10s data. From the differences to ATLAS we estimate the current accuracy at 50 mbar was 4-5%. A flight segment 
with data from Argus and ATLAS is shown in Figure 5. These result are in good agreement with a comparison of ER-2 and 
OMS balloon instruments during coincident measurements in Alaska (S. Alex, Harvard University, personal 
communication). 
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Figure 5: Time series of simultaneous measurements of N,0 of Argus and ATLAS aboard the ER-2 in November 1998. 

6. Future Modifications and Plans 

Our next field mission will be aboard the NASA ER-2 airplane during the NASA SAGE Ozone Loss and Validation 
Experiment (SOLVE). We plan to permanently adapt Argus for the ER-2 by changing the heat management to a forced air, 
convective cooling system. We will improve some key electronic components in the harmonic detection to reduce electronic 
noise. The laser temperature stability will be increased by switching from diode sensors to resistive elements and therefore 
reducing the temperature noise due to the AM detection by the diodes of high frequency signals radiated by the digital 
electronics. A major improvement of the precision down to 0.5% is expected by using a 36 m astigmatic Herriott cell, 
essentially doubling the path length Further we will improve the data analysis for the aircraft environment, and increase the 
data rate to 0.5 Hz. An extensive intercomparison of N 2 0 instruments, including ATLAS, in September 1999 will allow 
determination of the accuracy and precision of the modified Argus and assure the continuation of the ATLAS data series. 
With all these efforts we expect to reduce the accuracyof Argus to 2-3%. 


ATLAS/Argus comparison 981115: time series 
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of virtual walls, vve can product; functionality 
on microchips that is difficult to realize by 
other methods. For example, concentrating 
samples in microfluidic systems is nontrivial. 
Methods reported mclude sample stacking 
using “isoelectric focusing” (iff) and on-chip 
solid phase extraction ( 31 ). The virtual walls 
provide an intuitive method to concentrate 
liquids on microchips. Virtual walls can also 
mimic the function of lungs by exchanging 
components between liquid and gas phases. 
In combination with in situ-constructed stim- 
uli-responsive hydrogel components ( 25 , 32 ), 
more complex functions could be added to 
realize even greater control in microfluidic 
systems. We believe that the approaches we 
describe provide enhanced functionality and 
design flexibility that expand the toolbox for 
the design and fabrication of microfluidic 
systems. 


mation of SAMs in unv^anted areas. At the end of the 
process, flow of the siUne solution was always halted 
before stopping the how of pure solvent. The flow 
rates of solvent and solution were usually the same, 
either 1 or 2 ml/min; the flow time was 2 to 3 min. 
The channels were cleaned by sequentially flushing 
with 10 ml of hexane and 10 ml of methanol 
followed by drying with a stream of dean air or 
nitrogen. 

27. A pipette tip with inner diameters of 8 mm at the 
bottom and 6 mm at the top was fixed onto the 
channels. The maximum pressure that virtual walls 
can withstand was determined by gradually adding 
deionized water into the pipette and then measuring 
water height when a bulge developed. 

28. The surface free energy of Rhodamine B dilute solu- 


tion is lower than that of deionized water, re s ; 
a lower maximum pressure jj 
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The Detection of Large 
HN0 3 -Containing Particles in 
the Winter Arctic Stratosphere 
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Large particles containing nitric acid (HN0 3 ) were observed in the 1999/2000; 
Arctic winter stratosphere, These in situ observations were made over a targe- 
altitude range (16 to 21 kilometers) and horizontal extent (1800 kilometers), 
on several airborne sampling flights during a period of several weeks. With*: 
diameters of 10 to 20 micrometers, these sedimenting particles have significant' 
potential to denitrify the lower stratosphere. A microphysical model of nitric 1 
acid trihydrate particles is able to simulate the growth and sedimentation of. 
these large sizes in the lower stratosphere, but the nucleation process is not yet: 
known. Accurate modeling of the formation of these large particles is essential" 
for understanding Arctic denitrification and predicting future Arctic ozone- 
abundances. ^ 


'.4* 


Polar stratospheric cloud (PSC) particles play 
a well-recognized role in the chemical loss of 
stratospheric ozone over the Antarctic and 
Arctic regions in winter/early spring ( 1 - 3 ). 
PSCs are formed at low temperatures (<200 
K), characteristic of the winter stratosphere, 
through the co-condensation of water and 
HN0 3 . Heterogeneous reactions on PSC par- 
ticles produce active chlorine species. The 
sedimentation of PSC particles irreversibly 


each winter. NO y is the sum of pri| 
reactive nitrogen species, of which ^ 
NO, N0 2 , N 2 O v and C10N0 2 are impa 


removes HN0 3 (denitrification) and water 


(dehydration) from the lower stratosphere. 
Denitrification slows the return of chlonne to 
its inactive forms and, hence, enhances ozone 
destruction in a winter/spring season ( 4 - 6 ). 

Evidence of denitrification is abundant in 
both polar stratospheres from both m situ and 
remote observations ( 2 , 7 - 11 ). Denitrifica- 
tion is most intense over the Antarctic region, 
where large fractions of available NO y are 
irreversibly removed from the stratosphere 


in the lower stratosphere ( 12 ). In pnn£ 
denitrification must involve the sedi^ 
tion of large HN0 3 -containing particl^. 
14 ), but important details of the proces^ 
not been confirmed by observation o3 
descnbed theoretically. The concentranj 
denitrifying particles must generally be| 
less than that of the background liquid 3 
aerosols (CU-fam diameter). Too littl^ 
densable material (HN0 3 , H 2 0) exists|j 
stratosphere to allow all particles tow 
simultaneously to sizes sufficient forj 
sedimentation, particularly when ice doj 
form. 3 

Here we report observations in the ^ 
winter stratosphere of a newly identified^ 
of large HN0 3 -contaming particles. Th| 
were nh^prverl hprween JanuaTy 


£ 




1026 


9 FEBRUARY 2001 VOL 291 SCIENCE www.sciencemag.org 






particles contribute substantially to denitrifi- 
cation in the lower stratosphere? Third, how 
do these particles form in the lower strato- 
sphere? To address the first question, we used 
a statistical simulation o determine the size 
and number concentration of particles sam- 
pled in interval 2 (Fig. 2 ). This 800-s interval 
was chosen because the average NO signal 
due to particles is large (21 ppbv) and has a 
relatively constant standard deviation. This 
suggests that the size distribution is also rel- 
atively invariant over this interval. The large- 
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particle NO^ signal is defined as the differ- 
ence between the front and rear NO y values 
and thus only includes NO^ associated with 
particles larger than 2 pm (Fig. 2). The sta- 
tistical simulation accounts for the sampling 
of more than one particle in any 1 -s sampling 
period and for the spreading of HN0 3 from 
an individual particle over more than a single 


— ou 

NOy (ppbv) 

r ReSuttS from the statistical simulation of 
NC V difference values for interval 2 on 20 Janu- 
ary. The occurrence histogram of NO difference 
values is shown in (A) for the 800 s of observa- 
tions in interval 2 (solid line) and for 20,000 s in 
the simulation (dashed line). The larger simula- 
tion period reduces the variability in the simula- 
tion histogram in (A). The best -fit size distribution 
for interval 2 is given by the Gaussian functions N 
(# cm pm' 1 ) - 2.8 X 10’ 3 exp[-20(D - 
3.5) ] + 4.3 x 10 5 exp(-0.083(D - 14.5) 2 ] 
where 0 is the particle diameter in pun (Fig. 4) 
aussian functions are chosen for convenience to 
represent particle populations at large and small 
sizes. Lognormal distributions could also be used. 
The simulated time series is produced by super- 
imposing the released HNO, from all particle 
sizes and averaging it to 1-s'ntervals to match 
the sampling period of the NC) instrument. The 
differences between the simulated and observed 
histograms are shown in (B). The differences are 
very sensitive to the features of the large mode. 
For example, shifting the large mode by ±2 pirn 
would produce significantly larger differences in 
(BJ than shown. The small-size mode in Fig. 4 is 
required to match the histogram peak below 10 
ppbv The exact number and size of these smaller 
particles is more uncertain than for the larger 
particles. The small mode represents only a sub- 
set of particles present in the atmosphere below 
a 5 -tJ-m diameter. The histogram comparison is 
not sensitive to the presence of additional parti- 
cles with sizes between the two size modes and 
with concentrations CIO' 5 particles cm" 3 

pm 



10 15 20 

Particle diameter (pm) 

Fig. 4. Results from the statistical simulation of 
NO^ values in inteivai 2 on 20 January. Panels 
show (A) the adjusted best-fit size distribution 
and number concentration, (B) the equivalent 
gas-phase mixing ratio of HN0 3 derived from 
the distribution, and (C) the derived flux of 
HN0 3 at 60 hPa (—19.5 km). The size distribu- 
tion is given by the Gaussian functions in the 
caption of Fig. 3 adjusted to account for size- 
dependent sampling efficiency. The adjust- 
ments, which are in addition to a basic particle 
enhancement factor of 12.8 (72), vary from 
+ 30 to -10% over the size range and are, 
calculated with a fluid dynamical model of the 
particle separator. The concentration integral 
of the large (small) mode is 2.3 X 10~ 4 cm" 3 
(2 x 10 cm 3 ), with an estimated uncertain- 
ty of ±30%. The lower limit of the vertical axis 
in (A) corresponds to the detection of a single 
particle over the 800-s observation period The 
total gas-phase HNO a of the large (small) 
mode is 1.5 ppbv (0.2 ppbv). The HNO, flux of 
the combined distribution is 5 X 10 9 molecules 
cm km day' 1 or 2.2 ppbv km day' 1 at 60 
hPa. 
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1 -s period (20). The NO time > 
simulated by sampling a randomly <. 
population of NAT particles along ; 
fright track (Fig. 3). The numberY 
tions in each size interval are the 
parameters. The represen tativenes-: 
size distribution was determined b\ 
mg the occurrence histogram of N< 
ence values with that from the simui , 
senes (Fig. 3). The best-fit size di. r 
was found by iteratively adjusting t : 
bution to minimize the difference h 
(Fig. 3B). 

Two particle size modes, shown ;; 
are required in our analysis to repro : 
occurrence histogram of intern! 2. T 
mode has a mean diameter of 14.5 4 ' 
total number concentration of 2.3 . 
cm 3 . This mode produces most of ] 
values above 10 ppbv in Fig. 3. The J 
mode at 3.5 jam has a number conce ; 
of -2 X 10 -3 cm ~ 3 and is well si: 
from the larger mode by size bins wit ■ 
particles. For these concentrations, ; 
particle is sampled by the front inlet • 
imately every 1.8 s, and a small pa;: 
sampled every 0.2 s on average The 
mass contained in the large mode is 
lent to an ambient gas-phase value . 
ppbv along the flight track (Fig. 4 E , 
mass of the more numerous smaller n ■ 
— 0.2 ppbv. Because the simulation is i 
difference between the two NO nu • 
ments, the small mode in Fig. 4 repi . 1 
only the largest sizes in one or more p l 
modes present in the atmosphere be ( 
diameter of 5 pm (16). With lowd 
speeds, these smaller modes form sepr. 
from the large-particle mode and nearer 
sampling altitude, and may include pa: 
composed of STS. j 

To answer the second question, wit 
these particles contribute substantia 
stratospheric denitrification, we note tJ 
existence of large HN0 3 -containing pa 
means that denitrification has alread 
curred because particles sediment duri 
growth process. An instantaneous sedir 
tion flux of KN0 3 can be calculated.4 
simulated size distribution by convolvij 
fall speed and NO v content for each'pl 
size in 1-pm size intervals (Fig. 4C).|| 
speed of a 14-|pm spherical particle 
km day ’ 1 at 18 to 19 km (21). Thei 
flux for interval 2 is 5 X 10 9 moleculet; 
km day 1 ( 2.6 ppbv km day -1 ) aft 
altitudes. The flux values for manvS 
flight intervals on 20 January are siieuK 
the basis of average NO^ difference^ 
and occurrence histograms for thesejl 
vuis. Thus, some denitrification mu|g 
occurred above flight altitudes in thej® 
before 20 January in order to produgj 
large particles observed on that day if® 
widely distributed air parcels. The ^ 


1 


9 FEBRUARY 2001 VOL 291 RriFNrp 

butNCE www.sciencemag.org 




m 


cannot account for differs rices in observed 
particle number concentrations. Concentra- 
tions depend on the spatial and temporal vari- 
ation of the particle nucleation process 
throughout the low-tempt i ature regions of 
the vortex. 

The sensitivity of large particle formation 
to available HN0 3 and temperatures near ice 
saturation was examined with a one-dimen- 
sional model of aerosol growth and sedimen- 
tation (Fig. 6) {20). The remits show that the 
maximum NAT particle size that could be 
found at ER-2 cruise altitudes (50 to 70 hPa) 
is near 20 pirn in diameter tor average HN0 3 
values observed above ER-2 altitudes in mid- 
January (7 to 9 ppbv) (26). This maximum 
size is consistent with the NO^ observations 
(Fig. 4). At ER-2 sampling altitudes, a 25% 
reduction in available HN0 3 causes a 10% 
reduction in particle size (20 to 18 p-m). 
Similarly, a 2 K increase in temperature re- 
duces the particle diameter by <10%. Parti- 
cle sizes >22 p.m can only be produced with 
a uniform abundance of >10 ppbv HN0 3 
above the aircraft, which is greater than that 
observed in January to March 2000. These 
and larger panicles were therefore unlikely to 
be found at or above a 6()*hPa pressure alti- 
tude (—19 km). 

The particle growth model results show 
that denitrification by large particles is 
mostly controlled by the kinetics of HN0 3 
vapor deposition onto particles at condi- 
tions away from thermodynamic equilibri- 
um of the NAT condensed phase. This is a 
result of the low numbe concentrations of 
particles and their large fall speeds. The 
assumption that NAT particles are in equi- 
librium, as used in various ways in most 
current atmospheric model simulations of 
denitrification, will lead to incorrect simu- 
lations of denitrification Thus, the obser- 
vations and analysis presented here suggest 
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a need for fundamental changes in the vor- 
tex-scale modeling of denitrification. Rep- 
resentative modeling of denitrification will 
require accounting for particle growth 
along sedimentation trajectories, thereby 
convolving the temporal and spatial coor- 
dinates of the winter vortex in a detailed 
way. 

Critically lacking for representative model- 
ing of denitrification is an identified and veri- 
fied nucleation process for this newly identified 
class of large particles. If nucleation involves 
preexisting particles, the process must effect a 
high selectivity because background aerosol 
number concentrations are near 10 cm 3 in the 
lower stratosphere ( D > 0.008 p.m) (27, 28). 
whereas the computed large-particle concentra- 
tion Ln interval 2 is near 10" 4 cm -3 . Selectivity 
is the key to the growth of the large-particle 
population observed here and their denitrifica- 
tion potential. If all available background par- 
ticles grew and took up typical stratospheric 
HN0 3 amounts as NAT (or NAD), the average 
particle diameter would be <1 pun, and the 
sedimentation speed of individual particles 
would be <0.1 km day" 1 (/J). This speed is 
not adequate to explain significant denitrifica- 
tion as obsen/ed in both polar regions (13, 14, 
29). Selectivity could be the consequence of (i) 
a very low freezing rate of background particles 
to form NAT (or NAD) at temperatures below 
- 197 K (30. 31) t (ii) the unique composition of 
a few nuclei (32, 33), or (iii) ice formation (10, 
13, 34) possibly involving lee-wave clouds 

m 

Antarctic satellite observations show that 
denitrification occurs before dehydration as 
temperatures decrease in early winter (29, 35). 
The temporal separation of the HN0 3 and water 
removal processes could result if the large par- 
ticles reported here also form in the Antarctic 
stratosphere before the onset of ice saturation 
temperatures. The more extensive low temper- 


Fig. 6. Growth calculation esults 
for maximum NAT particle di- 
ameter as a function of sampling 
pressure, HN0 3 mixing ratio, and 
temperature (20). Particles are 
assumed to grow and sediment 
in a one-dimensional atmo- 
spheric model defined by a con- 
stant temperature of 188 or 190 
K ( an H 2 0 mixing ratio of 5 
ppmv, and a given HN0 3 mixing 
ratio (horizontal axis). Particle 
growth occurs wherever NAT su- 
persaturation conditions exist 
and, hence, begins at pressures 
between 13 and 23 hPa (25) for 
the highest and lowest HNO } 
mixing ratios, respectively, and 
extends down to the highest pressure shown (150 hPa). The solid curves indicate the constant 
HNO. mixing ratio and temperature values in the model that will produce the indicated particle 
diameter at a given pressure level (vertical axis). The horizontal dashed line is the ER-2 pressure 
level for sampling interva 2 in Fig. 2. NAD particles grow more slowly than NAT particles for the 
same HN0 3 partial pressure and, hence, would grow to smaller sizes (—20% less) for the same 
atmospheric conditions. 



atures in the Antarctic suggest that large 
cles are formed there every winter, * 
their formation in the Arctic may not oca 
year. Before the 1999/2000 Arctic observ 
dilute populations of large panicles m 
vonex may have eluded detection becaua 
able instruments did not sample the von 
the appropriate times and places. 

By the end of March 2000, large 
losses (up to 70% near 20 km) had oe; 
in the lower stratosphere throughout th j 
tic vortex (23. 36). The significant la 
denitrification observed throughout th; 
tex is expected to have enhanced the 
cal loss (id). Arctic ozone abundant, 
remain vulnerable to increased winter 
loss in the coming decades as anthrop : 
chlorine compounds are gradually re : 
from the atmosphere, particularly if 
concentrations of greenhouse gases i 
cooling in the polar vortex (37) and tri: 
increasing water vapor continue in the 
stratosphere (38). Both effects mcra - 
extent of PSC formation and, thereby . 
trification and the lifetime of active c 
(29). The role of large particles in efl 
denitrification in these future scena- 
likely quite important. Substantial am 
spread denitrification could possibly 
chemical ozone loss in the Arctic e / 
total chlorine levels fall below 2 p y\ 
2070 (10). In this case, recovery of : 
ozone would be delayed until chlorine 
decline further. Models of denitnficat: : 
accurately represent the formation o: 
particles will be required for accurate : 
tions and understanding of future ozc ' 
els in the variable Arctic vortex 
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00 With three instruments on board 
A ER-2 high-altitude ajrcraft as it 
Jhe lower stratosphere inside the Arc- 
el (Ftg- 1) (?*5, 16). Two instruments 
'.MO and HN0 3 , respectively, in 
gas phase and in aerosols, and one 
[fd the size and number concentration 
l^liiiua] particles. The presentation here 
^/primarily on the NO y instrument, 
ias two independent measurement 
hn order to distinguish sizes greater 
Indian ~2 p-m in diameter (“front” and 
jfesiO data, respectively) (Fig. 2). 
*fenost extensive population of large 
was observed at aircraft cruise alti- 
to 21 km) on 20 January 2000 
g^; flight between ICiruna, Sweden, and 
Jorth Poie (Figs. 1 and 2). Isolated peaks 
lont NO v data (Fig. 2C) at the begm- 
JjUd end of the flight represent the sam- 
'0qU individual HN0 3 -contaimng parti- 
throughout most of the flight, however, 
{writ NO^ data show a continuous senes 
$cs (Fig. 2D) due to the sampling of a 
eMjumber concentration of particles. The 
*^|uivalent sphere diameter associated 
^individual particle peak can be cal- 
B^'if typical PSC particle compositions 
burned (17). The most likely solid phas- 
ofm below 200 K are nitric acid tnhy- 
iHNOrsOHjO) (NAT) and nitnc acid 
3Jjte (HN0 3 *2H 2 0) (NAD) (3, 18, 19). 
gf^large particles are unlikely to be liquid 
j£rcWed ternary solution (STS) composed 
®)' 3 , H 2 0, and H 2 S0 4 (20). For a single 
m^-tomposed of NAT, the relation be- 
CtSkiiameter D and the recorded NO y val- 

K i(p.m) = 4.7 [NO, (ppbv)] 1 ' 3 (ppbv, 
_ t billion by volume). Thus, the data in 
Sg^suggest that particles of sizes much 
afe^than 2 |im in diameter were present 
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The NO y observations raise three ques- 
tions. First, how can the size and number 
concentration be determined for these large 
HN0 3 -containing particles? Second, do these 
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over large horizontal regions (1800 km) in 
the center of the Arctic vortex. The presence 
of large particles is confirmed by the direct 
HN0 3 and aerosol size measurements (16). 


Fig 1. Map of the ER-2 
ground track on the flight 
of 20 January 2000 (yel- 
low line) beginning in 
Kinjna, Sweden (67°N), 
reaching a maximum lat- 
itude of 89°N, and return- 
ing to Kiruna. The tem- 
perature field is shown by 
color shading and labeled 
contours (in K) on the 430 
K isentropic surface (—19 
km) at 1200 UTC on 20 
January. The edge of the 
Arctic vortex (red con- 
tour) is defined by a po- 
tential vorticity value of 
2.1 X 1G~ 5 K m 2 kg -1 
s‘ \ The thick portion of 
the yellow line indicates 
that portion of the flight 
track where large parti- 
cles were present. 


Fig. 2. Recorded data A 140 L ‘ ' ' ' ’ ' ‘ ' ‘ 20 January 2000 - 

from the ER-2 flight of ■ 2 I - 

mm 11 Hijhiiii i 

ues from the front z 40- i V^|H || jjllM IMLl I I H i 

(black) and rear (red) 20^ ^ 1 

inlets are shown for the p t . i ■ ■ ■ 2 

entire flight in (A), ^ A j\ *■ 20 

Temperature (black) ^ 205 A \ \ / \ / - 19 

and altitude (red) are * : l A/-^ \ / v \J .18 

shown in (B). The 800-s £ 200 - V 1 W 

irtervals labeled in (A) ~ • An * V i 1 

as "V (green) and ''2’' ^ 195- J Hr * 2 /'-A\ /V A Id ~ 16 

(blue) are expanded in £ : / V#Vl/M “ 15 

(C) and (D), respective- & J t . ■ , . 
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jjv-of the Arctic vortex in a given winter 
Sends on the spatial and temporal extent of 
Renting particle populations. If a flux 
jj]xf!>of 2.6 ppbv km day” 1 were to be 
gained for a period of only days near 20 
^jrf a column model, a significant fraction 
!$&) of available HN0 3 above 20 km 
jjtljd-be removed. 

j etwee n 20 January and 7 March, large 
^-containing particles were present in 
Bditional sampling flights inside the Arc- 
ortex. The number and spatial extent of 
cles observed on these later flights were 
|^ficant!y less than those observed tin 20 
ry. As the stratosphere warmed above 
isaruration temperatures by mid-March, 
particles no longer appeared in the NO 
|On these later flights, two sampling 
p contained particle populations that 
e|notably different from those observed 
j^ anuar ^‘ Periods occurred at alti- 
^etween 15 and 16 km. On 31 Janaary, 

f , burst of large particles was sampled 
m .n average NO y difference value many 
es^that of interval 2. With diameters esti- 
jcl to be near 20 jim, the instantaneous 
0^ flux far exceeds that found in intuwal 
' om the 3 February flight, the analysis of 
500-s interval shows a population of pri- 
y. 12-fj.m-diameter particles with con- 
iions near 10 3 cm -3 and an estim ued 
Comparable to that of interval 2. 
^nitrification was widespread in (he 
^2000 Arctic winter, as estimated with 
?| r measurements and the NO - 
correlation as a reference (7, 32). 

$3 f 1 cat ion exceeded 50% in many 
JCd air parcels in the 16- to 21 -km 
^altitude range throughout the vortex 
jBe; January to March period, including 
st:pf the 20 January flight. Widespread 
unification was also found in satel ite 
Rations (23). The unusually low . 
^pheric temperatures in this Arctic \ 
^beginning in December, increased < 
gel i hood of PSC formation and voru-x I 
Station (24), The different large- i 
“^'Populations described here all ccn- ° 
M to vortex denitrification through 
^formation and sedimentation. It is 
?£§* ve n the low winter temperatures, 
stinilar populations present before 20 
'2$ furTher contributed to observed Fi 
[location. 2 . 

Ir $ns wer to the third question, as to how 
grge particles form in the lower strato- di 
5- We note that NAT particles above 10 
^diameter found between 1 6 and 2 1 km du 
their growth and sedimentation nt dii 
^ al ly higher altitudes. The altitude frc 
^.ntent is a function of growth rate ant an 
action velocity. To demonstrate how S01 


Reports 

served in interval 2 (Fig. 5) Growing parti- 
cles are advected by atmospheric winds and, 
hence, will be transported horizontally and 
vertically while they gravitationally sedi- 
ment. The trajectory calculations take into 
account air parcel advection simultaneously 
with particle growth and sedimentation. The 
time evolution of particle size, pressure alti- 
tude, temperature, and NAT supersaturation 
ratio for particle sizes between 6 and 18 jim 
in diameter shows several important features 
(Fig. 5). First, the trajectories for all sizes 
terminate (reach zero size) in NAT-supersat- 
urated air, indicating that the observed sizes 
are consistent with atmospheric and micro- 
physical constraints of the trajectory model. 
Air masses encountered by particles must 
have temperatures that cause NAT supersat- 
uration (short times at higher temperatures 
can be survived if particles are large enough). 
Second, growth of the largest particle re- 
quires 6 days and an ~-6-krn vertical dis- 
placement at temperatures below NAT satu- 
ration values. Third, the trajectories terminate 
at different altitudes, temperatures, and NAT 
supersaturation values. The horizontal loca- 
tions of the termination points also vary (25). 
These preliminary results suggest that there 
are no obvious preferred locations or satura- 


tion conditions for nucleating the NAT parti- 
cles observed in interval 2. The accumulation 
of HN0 3 in large particles in interval 2 caus- 
es some denitrification over the range of ter- 
minal altitudes. Finally, the trajectory tem- 
peratures at the start of particle growth (188 
to 194 K) are above ice saturation by a few 
degrees. However, given the uncertainty in 
the meteorological analyses and atmospheric 
variability, ice saturation could have occurred 
along or near the trajectories and possibly 
played a role in particle nucleation. Particle 
trajectories also were calculated forward in 
time for interval 2. The results show that a 
20-p.m particle would continue to grow for 
another day, reach a maximum diameter of 22 
M>m at 150 hPa, and then evaporate (causing 
nitrification) within a few hours. Many in- 
stances of nitrification (NO, in excess of the 
reference value) were observed in the vortex 
on the ER-2 flights. 

Isolated particles sampled in interval 1 
were significantly smaller than the mean size 
found in interval 2 (Fig. 2). The results of a 
similar particle trajectory analysis for interval 
1 are consistent with the observed maximum 
size of —13 (xm in diameter. Although the 
trajectory analysis can account for such dif- 
ferences in size between different locations, it 
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.? panicle sizes in Fig. 4 evolved aftet 
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tu. nartirio in,. A- ^ , e ‘ and l™ HNO, supersaturation ratio with respect to NAT (5 1 at 

dSmete reiche were* TheT VerSUS - time (,n da ^' The trajectories terminate where the panicle 

ass s rc ssisussiir 




www sciencemag.org SCIENCE VOL 291 9 FEBRUARY 2001 




I 




wf'ciMS. air was sampled with an inlet specially 
oned to separate gas- and particle-phase HN0 3 
gaging a modified virtual impactor technique. Dur- 
most flights, gas and particle phases were sam- 
ecj -alternately for periods of -3 min each. The 
5 data confirm the NO y component of the large 
icies as HNOj on 20 January and provide size and 
concentration evidence consistent with the 
^^’observations. With the MA5P probe, the distri- 
fon of particles was measured for sues from 0.3 to 
M.m in diameter on 20 January. The probe data 
|, frn the presence of many small PSC particles 
i^een 0.3 and 2 p. m and a few larger particles up 
>0 p,m in diameter. The probe data for the size and 
nber of the larger particles are nominally consts- 
t* with the results in Fig. 4, although the statistical 
[ ehainty l5 high because of low count rates. The 
nber of particles in the size range of the small 
Ration mode in Fig. 4 is also nominally consistent 
h*the probe data. 

direct measurements of particle shape are avail- 
jh However, particles composed of solid HN0 3 
Irates are generally assumed to be nonsphencal. 
in at sizes of <S p,m, because of lidar depolariza- 
Vmeasurements [e.g .. (40)1- Although modest cor- 
tions to fall speed and sampling efficiency calcu- 
ons may ultimately be warranted because of non- 
lerical shapes, the assumption of sphericity is 
)6ted here throughout for simplicity, 
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The Indian Ocean Experiment: 
Widespread Air Pollution from 
South and Southeast Asia 

J. Lelieveld, 1 * P. J. Crutzen , 1,2 V. Ramanathan , 2 M. O. Andreae , 1 
C. A. M. Brenninkmeijer , 1 T. Campos , 3 C. R. Cass , 4 
R. R. Dickerson , 5 H. Fischer , 1 J. A. de Gouw , 6 A. Hansel , 7 
A. Jefferson , 8 D. Kley , 9 A. T. J. de Laat , 6 S. Lai , 10 
M. G. Lawrence , 1 J. M. Lobert , 2 O. L. Mayol-Bracero , 1 
A. P. Mitra , 11 T. Novakov , 12 S. J. Oltmans , 8 K. A. Prather , 13 
T. Reiner , 1 H. Rodhe , 14 H. A. Scheeren , 6 D. Sikka , 15 J. Williams 1 

The Indian Ocean Experiment (INDOEX) was an international, multiplatform 
field campaign to measure long-range transport of air pollution from South and 
Southeast Asia toward the Indian Ocean during the dry monsoon season in 
January to March 1999. Surprisingly high pollution levels were observed over 
the entire northern Indian Ocean toward the Intertropical Convergence Zone 
at about 6°S. We show that agricultural burning and especially biofuel use 
enhance carbon monoxide concentrations. Fossil fuel combustion and biomass 
burning cause a high aerosol loading. The growing pollution in this region gives 
rise to extensive air quality degradation with local, regional, and global impli- 
cations, including a reduction of the oxidizing power of the atmosphere. 


Until recently. North America and Europe 
dominated the use of fossil fuels, resulting in 
strong carbon dioxide emissions and global 
warming (7). The fossil energy-related C0 2 
release per capita in Asia is nearly an order of 
magnitude smaller than in North America and 
Europe (2). However, Asia is catching up. 
About half of the world’s population lives in 
South and East Asia, and hence the potential 
for growing pollutant emissions is large. In 
China, many pollution sources reduce air 
quality (3-5). In rural residential areas, nota- 
bly in India, the burning of biofuels, such as 
wood, dung, and agricultural waste, is a ma- 
jor source of pollutants (6). In urban areas, 
the increasing energy demand for industry 
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and transport propels fossil fuel utilization 

(7). 

Here we evaluate measurements of the 
Indian Ocean Experiment (INDOEX) to char- 
acterize the atmospheric chemical composi- 
tion of the outflow from South and Southeast 
Asia, from January to March 1999 during the 
dry winter monsoon (8). During this season, 
the northeasterly winds are persistent, and 
convection over the continental source re- 
gions is suppressed by large-scale subsi- 
dence, thus limiting upward dispersion of 
pollution (9). Our analysis is based on mea- 
surements from a C-130 and a Citation air- 
craft operated from the Maldives near 5°N, 
73°E, the research vessels Ronald H. Brown 
and Sagar Kanya , and the Kaashidhoo Cli- 
mate Observatory (KCO) on the Maldives 
(Fig. 1). During the campaign, the location of 
the Intertropical Convergence Zone (ITCZ) 
varied between the equator and 12°S. Hence, 
transport of primary pollutants and reaction 
products toward the ITCZ could be studied 
over an extended ocean area where pollutant 
emissions are otherwise minor. By perform- 
ing measurements across the ITCZ, the pol- 
luted air masses could be contrasted against 
comparatively clean air over the southern 
Indian Ocean. Furthermore, we used the mea- 
surements to evaluate the numerical represen- 
tation of these processes in a chemistry gen- 
eral circulation model (GCM) (10). The mod- 
el was subsequently applied to calculate the 
large-scale atmospheric chemical effects ot 
the measured pollution. 

Aerosol chemical and optical measure- 
ments were performed from both aircraft, the 
R/V Brown , and KCO. The latter is located 
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Severe and Extensive Denitrification in the 1999-2000 
Arctic Winter Stratosphere 


P.J. Popp 1,2 , M.J, Northway 1 ’", J.C. Holecek 1 ’ 2 , R.S. Gao 1 , D.W. Fahey 1,2 , J.W. Elkins 3 , 
D.F. Hurst 3,2 , P.A. Romashkin 3 ' 2 , G.C. Toon 4 , B. Sen 4 , S.M. Schauffler 5 , RJ. 
Salawitch 4 , C.R. Webster 4 , R.L Herman 4 , H. Jost 6 , T.P. Bui 6 , P.A. Newman 7 , L.R. Lait 7 


Abstract. Observations in the 1999-2000 Arctic winter strato- 
sphere show the most severe and extensive denitnfication ever 
observed in the northern hemisphere. Denitrification was inferred 
from in situ measurements conducted during high-altitude aircraft 
flights between January and March 2000. Average removal of 
more than 60% of the reactive nitrogen reservoir (NO y ) was ob- 
served in air masses throughout the core of the Arctic vortex. 
Denitrification was observed at alti.udes between 16 and 21 km, 
with the most severe denitnfication observed at 20 to 21 km. Ni- 
trified air masses were also observed, primarily at lower altitudes. 
These results show that denitrification in the Arctic lower strato- 
sphere can approach the severity and extent of that previously ob- 
served only in the Antarctic. 


Introduction 

Denitrification, defined as the permanent removal of NO y 
from an air mass by the sedimentation of polar stratospheric 
cloud (PSC) particles, can slow chlorine deactivation and thereby 
enhance polar ozone loss [Fahey ei al. , 1990; Rex et al., 1997]. 
Polar stratospheric clouds contain nitric acid (HNO}), the princi- 
pal component of stratospheric NO y . Severe (>50% NO y re- 
moval) and geographically extensive denitrification in the Ant- 
arctic lower stratosphere is a well-known phenomenon [Fahey et 
al, 1990; Fahey et al., 1989a; Santee et al., 1999]. The degree to 
which the Arctic stratosphere denitrifies is more variable, how- 
ever, due to generally warmer temperatures that result in shorter 
PSC lifetimes [Fahey et al., 1990; Santee et al., 1999; Tabazadeh 
et al., 2000]. Thus far, observations of denitrification in the Arc- 
tic vortex have been limited to those made by satellites and infre- 
quent in situ sampling from aircraft and balloon platforms [Fahey 
et al., 1990; Santee et al., 1999; Sugtta et al., 1998; Kondo et al., 
2000]. While severe denitnfication has been inferred from the in 
situ measurements for localized regions of the Arctic, [Fahey et 
al., 1990; Sugita et al., 1998], these observations are not consid- 
ered representative of the entire vortex. Although NO y removal 
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that is both severe and widespread has not been observed previ- 
ously in the Arctic, the region is arguably at the threshold for 
such denitrification [Waibel et al., 1999]. The observations pre- 
sented here show the 1999-2000 Arctic winter crossing the 
threshold for Antarctic-like denitrification. 

Results and Discussion 

In situ measurements of NO y , by catalytic reduction/chemi- 
luminescence [Fahey et al, 1989b], and nitrous oxide (N 2 0), by 
gas chromatography [Romashkin et al, 2001], were conducted in 
the Arctic vortex during 1 1 flights of the NASA ER-2 aircraft 
between 20 January and 12 March 2000. These flights were per- 
formed as part of the joint SAGE III Ozone Loss and Validation 
Experiment (SOLVE)/Third European Stratospheric Experiment 
on Ozone (THESEO) 2000 campaigns in Kiruna, Sweden (68° N, 
20° E). Denitrification is quantified using the correlation between 
the measured NO y and N 2 0 mixing ratios (reported every 70s) as 
shown in Fig. la. The reference value of NO y , NO y *, is given by 
the NO y -N 2 0 correlation established using remote measurements 
obtained by the MkIV balloon-borne interferometer [Toon, 1991] 
in the 1999-2000 Arctic vortex prior to the onset of denitrification 
(solid line in Fig. la). 

The aircraft data show severe denitrification between late 
January and mid-March, at altitudes from 17 to 21 km, and for 
N 2 0 levels between 40 and 170 ppbv (Fig. la). All air masses 
observed with N 2 0 levels less than 125 ppbv (encountered on 10 
of 11 ER-2 flights) were denitrified. The most severely denitri- 
fied air masses showed residual NO y levels (NO y remaining in a 
denitrified air mass) of less than 3 ppbv, resulting from more than 
80% of the NOy being removed. Data points above the NO y * ref- 
erence line in Fig. la indicate nitrification, thought to be caused 
by the evaporation of PSC particles that have sedimented from 
higher altitudes. To exclude measurements made when NO y - 
containing PSC particles were present in the sampled air mass, 
the in situ data shown in Fig. la (and throughout) represent ob- 
servations made at temperatures greater than the threshold tem- 
perature for the formation of nitric acid trihydrate (T^ a j). Val- 
ues of T nat were calculated continuously using in situ measure- 
ments of NO y (as equivalent HNO 3 ), water vapor, and pressure. 
Excluding these measurements ensures that NO y is irreversibly 
removed from the denitrified air mass, and not partially seques- 
tered in particle form. Residual NO y levels as low as 0.9 ppbv 
were also observed, although these measurements were made in 
air masses below PSC formation temperatures (data not shown). 

Remote measurements of NO y reported for a flight of the 
MkIV interferometer on 15 March 2000 show up to 40% denitn- 
fication at 19-20 km, in comparison to the 3 December 1999 
flight used to establish the NO y reference state (Fig. lb). These 
denitrifications are less than those measured onboard the ER-2, 
particularly for values of N 2 0 < 100 ppbv. Although the 15 
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Figure 1. (a) NOy-^O scatter plot for data measured during 1 1 
flights of the NASA ER-2. Data are colored as a function of the 
altitude at which the individual measurements were made. Un- 
filled squares represent measurements made by the MkIV bal- 
loon-borne interferometer prior to the onset of denitrification, 
flown inside the Arctic vortex from Esrange, Sweden (68° N, 20° 
E) on 3 December 1999. The NO y * reference state was deter- 
mined from a third-order polynomial fit to the MkIV data, given 
by 

NO y * = 17.3 - (0.0222-N 2 O) - (9.85- 10' 5 ) (N 2 O) 2 

for N 2 O values between 33 and 321 ppbv The dotted line repre- 
sents 50% of the value of NO y *» indicating the difference between 
moderate and severe denitrification, (b) NOy-^O scatter plot 
for data reported by the MkIV interferometer during balloon 
flights on 3 December 1999 and 15 March 2000. Data are col- 
ored as a function of the altitude at which the measurements were 
reported. The black line represents the NOy\ reference state 
shown in panel (a), (c) Composite NOy-^O scatter plot for the 
ER-2 data shown in panel (a) (blue symbols). Similar treatments 
of the data are shown for 7 flights of the Airborne Antarctic 
Ozone Experiment in 1987 (red symbols) and 8 flights of the Air- 
borne Arctic Stratospheric Expedition in 1988-89 (green sym- 
bols). Bins showing severe denitrification are represented by un- 
filled symbols. Bars show the standard deviation of the data 
within each N 2 O bin, and are not intended to represent the un- 
certainty in the mean NO y vdue. 


March MkIV 2000 observations were obtained at high equivalent 
latitudes, the tracer-potential temperature relationships observed 
during this flight were more representative of the vortex edge 
than the vortex core, which might explain the lower values of de- 
nitrification inferred from the MkIV measurements. Alterna- 
tively, the differences could arise from the different sampling 
characteristics (i.e. vertical resolution, altitude ceiling) of the re- 
mote and in situ data. 

Mixing can produce air masses with NOy-^O relationships 
that deviate from the NOy 1 ' reference state because of the non- 
linearity of the NOyN 2 0 correlation at low values of N 2 0 [Rex 
el al., 1999]. These air masses can result from isentropic mixing 
of air across the edge of the vortex, or from isentropic mixing of 
air masses within the vortex that have undergone non-uniform de- 
scent. Without knowledge of such mixing, air parcels which de- 
viate from the reference state could be described incorrectly as 
having been denitrified. The initial relationship between NO y and 
N 2 O measured inside the Arctic vortex by the MkIV interfer- 
ometer on 3 December. 1999 (Fig. lb) shows that the near-linear 
portion of the inilial reference state in the Arctic vortex extends to 
lower values of N 2 O (-20 ppbv) than had previously been as- 



Equivalent Latitude (deg. from equator) 

Figure 2. (a) Denitrification as a function of equivalent latitude 
observed during 1 1 flights of the NASA ER-2. Positive values of 
denitrification represent NO y removal. Unfilled black symbols 
represent mean denitrification values for measurements made at 
altitudes between 18 and 21 km. Error bars show the standard 
deviation of the data within each equivalent latitude bin. The 
dashed line at 65.2° N represents the mean equivalent latitude of 
the vortex edge (defined by potential vorticity) during the 1 i ER- 
2 flights, (b) Same as Fig. 2a for 6 Antarctic vortex flights of the 
NASA ER-2 between August 30 and September 22, 1987. Un- 
filled black symbols represent mean denitrification values for 
measurements made at altitudes between 16 and 21 km. The 
dashed line at 64.9° S represents the mean equivalent latitude of 
the vortex edge dunng the 6 ER-2 flights. 
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Figure 3. Vertical profile of redistributed NO y in the 1999-2000 
Arctic vortex. Positive values of denitrification represent NO y 
removal. Blue and red symbols represent measurements made at 
greater than 70° N and less than 70° N equivalent latitude, re- 
spectively. The solid black line representing NO/ was deter- 
mined from altitude bin-averaged N 2 0 values for measurements 
made throughout the winter at equivalent latitudes greater than 
70° N, This term represents the maximum available NO y at a 
given altitude. Maximum ER-2 flight altitudes of 20-21 km rep- 
resent a potential temperature of approximately 460-475 K. 

sumed (-40-50 ppbv) by Rex et ai 1 1999]. This feature makes it 
unlikely that transport and mixing substantially altered the NO/ 
reference state from that shown in Figs, la and lb. The in situ 
data set of long-lived trace gases obiained by the Observations of 
the Middle Stratosphere (OMS) balloon gondola during two 
flights above Kiruna also indicate that mixing did not signifi- 
cantly alter the NO y -N 2 0 relationship during this winter. 

Studies of PSC climatology show that consistently cold tem- 
peratures characteristic of the Antarctic vortex are highly favor- 
able for the formation and sedimentation of PSCs [Poole and 
Pitts , 1994]. Temperatures in the Arctic winter stratosphere, 
however, are typically wanner and more variable from year to 
year, resulting in a lower probability of PSC formation and deni- 
trification [Poole and Pitts, 1994]. These temperature patterns 
are consistent with previous in sit i measurements of NO y and 
N 2 0 [Loewenstein et ai, 1989] that show more extensive denitri- 
fication for the 1987 Antarctic winter than for the 1988-89 Arctic 
winter (Fig. 1c). While the 1988-89 Arctic data show evidence of 
only moderate NO y loss, the 1987 Antarctic data set indicates the 
onset of severe denitrification at N >0 levels less than 130 ppbv. 
Satellite data collected by the Microwave Limb Sounder (MLS) 
in more recent years also show more denitrification in the Ant- 
arctic vortex compared to the Arctic [Santee et ai, 1999]. 

When the 1988-89 and 1999-2000 Arctic data are compared, 
there is a significant difference in the inferred denitrification (Fig. 
Ic). In 1999-2000, denitrification occurs over a broader range of 
N 2 0 values, the onset of NO y loss occurs at higher N 2 0 values 
(190 ppbv), and mean residual NO, levels fall below 5 ppbv at 
the lowest values of N 2 0. In comparison to the Antarctic data 
set, the 1999-2000 Arctic data reveal denitrification over a 
broader range of N 2 0 values, and slightly more loss of NO y be- 
tween N 2 0 values of 1 10 and 170 ppbv. Only at the lowest val- 


ues of N 2 0 common to all three data sets (90-1 10 ppbv) does the 
Antarctic data show more removal of NO y . Thus, the severity of 
denitrification in the 1999-2000 Arctic vortex is comparable to 
that inferred from previous in situ measurements in the Antarctic 
polar vortex. The high values of NO y between 230 and 290 ppbv 
N 2 0 in the 1999-2000 Arctic data are due to the nitrified air 
masses at altitudes between 15 and 17 km (see Fig. la). 

In addition to causing denitrification, Antarctic temperatures 
are consistently low enough to cause dehydration via the sedi- 
mentation of ice particles. Satellite observations in the Antarctic 
vortex suggest that denitrification occurs before dehydration as 
temperatures decrease during the early winter [ Santee et ai, 
1995] In situ measurements of water vapor [May, 1998] were 
obtained during all flights of the NASA ER-2 in the 1999-2000 
Arctic vortex. Denitrification was not accompanied by any sig- 
nificant dehydration in this data set except briefly during the 
flight of 27 January 2000, consistent with the permanent removal 
of NO y at temperatures above the formation threshold for ice 
particles. 

The severe denitrification in the 1999-2000 Arctic stratosphere 
was observed throughout the core of the vortex. In Fig. 2a frac- 
tional denitrification is shown as a function of the equivalent 
(vortex-normalized) latitude [Nash et ai, 1996] of the measure- 
ments An equivalent latitude of 90° corresponds to the center of 
the vortex. Mean denitrification values, at altitudes between 18 
and 21 km, indicate that more than 60% of the expected NOy was 
removed from air masses sampled at equivalent latitudes between 
70° and 85° N (black symbols in Fig. 2a). The most severely de- 
nitrified air masses, with more than 80% of the expected NO y re- 
moved, were observed poleward of 73° N. Negative values in 
Fig. 2a indicate nitrified air masses, observed over the same range 
of equivalent latitudes but at lower altitudes. A similar treatment 
of the data for in situ measurements obtained in the 1987 Antarc- 
tic vortex is shown in Fig. 2b. The most severely denitrified air 
masses observed in the Antarctic show more than 90% of the 
NO y being removed at equivalent latitudes greater than 65° S, 
with mean denitrification values between 16 and 21 km (black 
symbols in Fig. 2b) exceeding 80% poleward of 75° S. Figure 2 
shows a striking similarity between the severity and extent of de- 
nitrification in the Arctic and Antarctic vortices in the observed 
years, with the Antarctic data set showing more severe denitrifi- 
cation over a greater range of equivalent latitudes. The apparent 
denitrification (and nitrification) at equivalent latitudes less than 
55° S in Fig. 2b occurs in air masses containing low values of 
NO y , for which small deviations from the NO/ reference state 
can cause large fractional denitrification. 

Severe denitrification over broad regions of the 1999-2000 
Arctic vortex was not detected in MLS satellite measurements 
[Santee et ai, 2000]. The MLS data suggest that only 20% of the 
available HNO3 was removed between equivalent latitudes of 70° 
and 80° N. Since the HN03/NO y ratio generally exceeds 80%, 
percent denitrification is approximately equal to that of percent 
HNO3 removal. The difference in the reported severity of denitri- 
fication in the 1999-2000 Arctic vortex may be due to the loca- 
tion, timing and vertical resolution (6 km) of the MLS HNO3 
measurements [Santee et ai, 2000], which causes denitrification 
to be reported as an average of many air parcels over a range of 
altitudes. 

A vertical profile of denitrification in the 1999-2000 Arctic 
stratosphere reveals that NO y removal was most intense in the 
vortex core at altitudes of 20-21 km (blue symbols in Fig. 3). Se- 
vere denitrification may have also occurred above 21 km, the 
maximum altitude attained by the ER-2. Below 21 km, maximum 
NOy removal decreases with altitude down to 16 km, below 
which no denitrification was observed. This altitude dependence 
results from lower amounts of available NO y at lower altitudes, as 
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indicated by the black line denoting NO y * in Fig. 3, and warmer 
temperatures (and thus decreased PSC activity) at lower altitudes. 
In addition, the denitrifica! ion observed in any given air mass re- 
flects a balance between the NO y removed from the air mass and 
any nitrification that has occurred. This conclusion is supported 
by the observation of some air masses between 15 and 17 km 
with higher than expected NO y mixing ratios. Such nitrification 
has been observed previously in the Arctic polar vortex by both in 
situ [ Hubler et ai , 19901 and remote methods [Kondo et ai, 
2000]. 

Implications 

Modeling studies have demonstrated that denitrification can 
significantly enhance ozone loss in polar spring [Satawitch et ai, 
1993; Rex et ai, 1997; Waibel et al., 1999; Tabazadeh el ai, 
2000]. The enhancement s expected to be greatest in denitrified 
air masses for which temperatures are too high to maintain active 
chlorine levels via heterogeneous reprocessing [Portmann et al., 
1996]. Consequently, denitrification may play a more prominent 
role in ozone loss in the warmer Arctic spring, despite the typi- 
cally higher levels of obseived denitrification reported in the Ant- 
arctic [Portmann et al., 1996]. Severe ozone loss was observed 
in the 1999-2000 Arctic winter stratosphere [Sinnhuber et al., 
2000; Richard et al. , 2001 ]. Modeling studies have partly attrib- 
uted this enhanced ozone loss to a delay in chlorine deactivation 
caused by denitrification [Sinnhuber et ai, 2000]. The 1999- 
2000 Arctic vortex experienced temperatures below the formation 
threshold for PSC panicles over a greater area than in any previ- 
ously observed winter [Manney and Sabutis , 2000]. These cold 
temperatures were accompanied by denitrification in the lower 
stratosphere approaching the severity and extent of that previ- 
ously observed only in the Antarctic. If the observed trend to- 
wards a colder Arctic stratosphere continues [Shindell et ai, 
1998], severe and extensive denitrification and the associated en- 
hanced ozone loss is likely. The observations in the 1999-2000 
Arctic vortex provide an important example of the denitrification 
that may become more frequent in future colder Arctic winters. 
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